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Abstract 

We present a two-Higgs-doublet model, with a Z3 symmetry, in 
which CP violation originates solely in a soft (dimension-2) coupling 
in the scalar potential, and reveals itself solely in the CKM (quark 
mixing) matrix. In particular, in the mass basis the Yukawa inter- 
actions of the neutral scalars are all real. The model has only eleven 
parameters to fit the six quark masses and the four independent CKM- 
matrix observables. We find regions of parameter space in which the 
flavour-changing neutral couplings are so suppressed that they allow 
the scalars to be no heavier than a few hundred GeV. 
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1 Introduction and notation 



One of the conceptually simplest extensions of the Standard Model (SM) 
of the electroweak interactions consists in allowing for uh > 1 gauge-SU(2) 
"Higgs" doublets. In such multi-Higgs-doublet models (MHDMs) CP viola- 
tion may occur in various places: in the quark mixing matrix (CKM matrix) 
just as in the SM, in the Yukawa couplings of the scalars to the quarksp 
in the mixing of the scalars (in particular, scalar-pseudoscalar mixing), or 
in the self-interactions (cubic and quartic interactions) among the scalars. 
Unfortunately, MHDMs in general lead to the existence of flavour-changing 
neutral currents (FCNC)I which are severely restricted by the experimental 
data. 

The simplest MHDMs are, of course, two-Higgs-doublet models (2HDMs) 
[1], which have lately been the object of intense scrutiny [2j. The Yukawa 
interactions of the quarks in the 2HDM are written 

2 

Cvnk = -Ql ^ (^(pk^kriR + (pkAkPRj + H.c, (1) 
fc=i 

where 0i,2 are the two scalar gauge-SU(2) doublets, (pk = iT2(pl for k = 1,2, 
Ffc and are (in general, complex) 3x3 matrices in flavour space, and 
Ql, nji, and pji denote the 3- vectors (in flavour space) of quark left-handed 
doublets, right-handed charge —1/3 quarks, and right-handed charge -1-2/3 
quarks, respectively. In order for the U(l) gauge group of electromagnetism 
to be preserved, the Higgs doublets are assumed to have vacuum expectation 
values (VEVs) of the form 

(0 10,1 0) = ( ^, ° ), (0 4>k0)=( ^'^'q ) , (2) 



with real and non-negative Vk- The quark mass matrices are then 
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= Y^Vke^'Tk, (3) 

k=l 
2 

Mp = Y.^kC^'^Ak. (4) 



fc=i 



-'^In this paper we neglect the lepton sector. 

^More precisely, quark-flavour-changing Yukawa interactions of the neutral scalars. 
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These are bi-diagonalized as usual by unitary matrices 

Ul^MJJl = Md = diag (m^, m„ m^) , (5) 
Ul^MpUl = M„ = diag (m„, m^, m^) , (6) 

and the CKM matrix is ^ = U^^U^- The quantity v = ^/vf + v"^ 
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{2^/2Gf) ~ 174 GeV is responsible for the masses of the and Z° 
gauge bosons. It is convenient to use the 'Higgs basis', 



= {vie-''^(Pi + V2e-''^(t)2)/v (7) 
G+ 

H2 = {v2e-''^ct>r-vre~''^ct>2)/v (9) 
{H + iA) 1^2 



(8) 



(10) 

in which only Hi has VEV, which is precisely v. The fields and G^ are 
the would-be Goldstone bosons. The field G^ is a physical charged scalar. 
The neutral fields /i, if, and A in general mix to form the three physical 
neutral scalars of the 2HDM. We define the matrices 

iV„ = v^e^'^Ti-vie^'-^T^, (11) 

Np = V2e-'''Ai-vie~'''A2, (12) 

and 

= U^^Nr.U'^, (13) 

AT. = Ul^N.Ul (14) 

Equation ([T]) then becomes 

-^Yuk = -diMddR - ulMuUr 

h H 

7k' {d-LMddR + ulMuUb) j=- [dLNddn + ulNuUr) 

iG^ - iA - 

-=- [diMadR - ulMuUr) -=- (dLNddn - ulNuUr) 

V2v v2f 

+ — u {M^V^L - VMdiR) d + —u {NlV^L - VNdiR) d 

V V 

+H.C., (15) 
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where d and u denote the column vectors in flavour space of the charge 
—1/3 and charge +2/3 quarks, respectively, in the mass basis, and 'yL,R are 
the chirality projection matrices in Dirac space. Since the matrices and 
Nu are not necessarily diagonal, the terms diNddR and ulNuUr in general 
include potentially problematic FCNC. 

We spot in the Yukawa interactions of equation f lT5]) three possible man- 
ifestations of CP violation: 

The CKM matrix V may contain a complex phase, just as in the SM. 

The matrices Nd and A^^ may be complex. 

The scalars h and H may mix with the pseudoscalar 

One further manifestation of CP violation may occur in the cubic and quartic 
interactions among the scalars. It is the purpose of this paper to present 
a 2HDM with an additional symmetry such that only the first one of the 
above four manifestations of CP violation occurs; namely, the matrices Nd 
and Nu are real, the scalars do not mix with the pseudoscalar, and the cubic 
and quartic interactions among the (neutral and charged) scalars respect CP 
invariance. Additionally, our model shows that the FCNC may be quite 
suppressed even when all the scalars have relatively low (less than 1 TeV) 
masses. 

2 The model: Yukawa couplings 

Our model is a 2HDM supplemented by a particular Z3 symmetry and by the 
usual CP symmetry. Let uj = exp(2i7r/3). Then, under the Z3 symmetry, 
the following matter fields transform as 

Qli uj'^Qli, Ql2 ujQl2, (16) 

nR3 UnR3, Pr3 UPR3, 

''if this mixing exists, i.e. if the three physical neutral scalars are mixtures of all three 
h, H , and A, then A is not a physical particle and it does not make sense to separate the 
physical neutral scalars into two scalars and one pseudoscalar. 
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and all other fields remain invariant. This symmetry forces the Yukawa- 
coupling matrices to have the following form |3]: 



Fi, Ai ~ X , ~ , A 






where the symbol x denotes a non-zero matrix entry. The standard CP 
symmetry forces all those non-zero entries of the Yukawa-coupling matrices 
to be real. Therefore, the mass matrices end up being 

e'^x \ 
M„ = e^^i ( a , (17) 
h c e^^y J 

e~'^a' \ 
Mp = e"*^^ I e-'%' e-'^c' x' , (18) 
y' J 

where 9 = 02 — Oi and a, b, c, x, and y' are real. In equation f[T7|) we 
have already assumed a rotation between nm and nij2 which renders zero 
the (1,2) entry of in the same way, in equation (fT8|) a rotation between 
Pfii and pji2 has been used to make (Ai)32 = 0. The matrices parameterizing 
the Yukawa couplings of H2 are 

-e^'^x/r 
Nn = e'^' I ra I , (19) 

rb rc —e^^y/r 



Np = e-'"' I -e''%'/r -e-'^c'/r rx' | , (20) 
ry' 

where r = V2/V1. 
Let 





Mrf, (21) 
M„, (22) 
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where O^'ji are real orthogonal matrices. It is then clear that 

Ui = e-diag(^e-|^e-«|^e-)0S. (24) 



\ \a'x'\ \a'b'x'y'\ J 

The CKM matrix is 

y = Of diag(l,e-,±e-)02, (27) 

where 

^ a'x'aybx ^^^^ ^ a'x'b'y'bx ^.3^^ 

I a'x'ayfca; I ' |a',T'6'?/'6,T| 

One sees that the complexity of the CKM matrix originates exclusively from 
the phase 39, which is the only phase with physical consequences in our 
model. One easily finds the matrices parametrizing the non-diagonal Yukawa 
couplings: 

^ Of { r\a\ I O^, (29) 



— \x\/r 





n ] 










r\a\ 




r\b\ 


r c| 


-\y\/r ) 










— \a'\/r \ 


-\b'\/r 


-\c'l 


jr r\x' 


) 


i^\y'\ 











Nu = of I -\b'\/r -\c'\/r r\x'\ I O^. (30) 



These matrices are real. Thus, in our model there is no CP violation from 
the FCNC matrices. 
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3 The model: scalar potential 



The scalar potential of our model is 

V = V,ym + VsB, (31) 

+ Y + Y (<^202)' + A3 0^02 + A4 0l02 0^01, (32) 

^SB = -|/i3|fe-^Vi02 + e^>^0i), (33) 



where V^ym respects the ZIs and CP symmetries of the model while Vsb breaks 
both those symmetries, but only softly. The soft-breaking term is unique and 
is as general as possible. Note that Vgym coincides with the Peccei-Quinn 
potential [4]. The minimization of the potential leads to the vacuum phase 
6 being equal to the phase in Vsb- Thus, in our model the origin of CP 
violation lies exclusively in a soft term in the scalar potential. 
The equations for vacuum stability read, besides 9 = {}, 

t^i = lyUsI — - Aii;i - (A3 + A4) 1^2, (34) 
ti'2 = l/^sl — - A2t^2 - (A3 + A4) vl- (35) 

V2 



If we define 



1/^3 I 



= , m^^ = rUj^ — X^v , (36) 



^ Note that ours is a model with soft CP breaking — the Lagrangian does not enjoy 
CP symmetry because of the presence of the /is term. This is distinct from a model [5] 
in which spontaneous CP violation is achieved through the addition to the Lagrangian 
of a soft (dimension-2) term which breaks some other internal symmetry but does not 
break CP. Spontaneous CP violation usually leads to CP violation in the scalar sector, in 
particular through scalar-pseudoscalar mixing. However, recently a model was found [6] 
in which there is spontaneous CP violation but the scalar sector still preserves CP. 
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then we easily find that the part of V which is bihnear in the fields is 

bilinear = ^ {A' + H') + m'^C- 

^ Xivf + X2VI + 2 (A3 + A4) vfyj ^2 

+ [Ai + A2-2(A3 + A4)]^i/' 

^^y^y^^^-!->^2vl + {Xs + X,)ivl-v!) ^3^^ 

V 

One sees that A does not mix with h and H. In our model there is no 
scalar-pseudo scalar mixing. 

Moreover, in our model there is no CP violation in the self-interactions of 
the scalars. This follows from the fact that in a general 2HDM there is only 
one gauge-invariant vacuum phase — 9 — and in our specific 2HDM there are 
only two terms in the scalar potential — those with coefficient \^^ \ exp {-^id) — 
which are sensitive to that phase. The vacuum phase adjusts in such a way 
as to offset the phase of those terms in the scalar potential so that the final 
potential has no phase at all. 



4 The fit: procedure 
4.1 First stage 

As seen in equations ( 12T|) and ( 122|) . the six quark masses depend only on ten 
parameters: |a|, |c|, ||/|, |a'|, |c'|, and \y'\. Then, from equa- 
tion fl27p . the CKM matrix V , which contains four independent observables, 
depends on one additional parameter, the phase 9% One thus has to fit ten 
observables by means of eleven parametersjj 

We have assumed throughout that the contributions to quark decays from 
tree-level diagrams with intermediate scalars are much smaller than the con- 
tributions from diagrams with intermediate W^. We thus assume that the 

^The CKM matrix additionally depends on the signs of a'x'bxay and of a'x'bxb'y', as 
seen in equation ([28| . 

^Even when the number of parameters is larger than the number of observables to be 
fitted, obtaining a good fit is not always possible. The fact that our model passes this test 
is interesting in itself. 
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SM extractions of \Vus\, \Vcb\, and \Vub\ still hold in our model. These three 
CKM-matrix elements and the quark masses are allowed to take any value 
within their Particle Data Group (PDG) allowed ranges In our fits \Vtd\ 
is left free, but we have found that, once the various experimental constraints 
to be discussed below are included, a good fit is obtained only when \ Vtd\ lies 
roughly in the SM-allowed range. 

We then proceed to analyze the FCNC of our model. These are governed 
by the matrices and Nu in equations ( 129|) and (!30|) . respectively. Those 
matrices involve the extra parameter r = V2/V1. 

In our analysis of the FCNC, we consider only their contributions to the 
mixing in the neutral-meson-antimeson systems K, B^,, Bg, and D. The rel- 
evant quantity is the off-diagonal matrix element M12 connecting each meson 
to the corresponding antimeson. That matrix element receives contributions 
both from an SM box diagram and a tree-level diagram involving the FCNC. 
We denote the latter by NP (for "New Physics" ) and write 

M12 = Mff + Mf/. (38) 

In order to shorten our text we shall follow the notation in the text- 
book |8] and freely use its equations with the prefix BLS. For the K system, 
Mf^ and the quantities relevant for its determination can be found in equa- 
tions (BLS-17.14), (BLS-17.16), (BLS-B.15), (BLS-B.16), and (BLS-13.50); 
the expressions for the other neutral-meson systems are obtained by straight- 
forward modifications of the quarks and mesons involved. The quark masses 
and CKM-matrix elements utilized in the calculation of M12 are those pro- 
duced by each of our fits; in addition, we use some other quantities shown in 
Appendix A. 

The calculation of is in equation (BLS-22.76)0 This calculation 

requires equations (BLS-22.29), (BLS-22.33), and (BLS-22.73). Since there 
is no scalar-pseudoscalar mixing in our model, one has 

Mr = <2 + Mi^', (39) 

where 

originates in the tree-level exchange of the pseudoscalar (parity-odd) 
A; 

^That equation contains a sign mistake in the hadronic matrix elements in the vacuum- 
insertion approximation, which we have corrected. 
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originates in the tree-level exchange of the two physical parity-even 
scalars 5*1 and S2, with masses mi and 1712, respectively. 

The scalars are mixtures of H and h through 



H 

h 



cos ip sin ip 
— sin ip cos ip 



Si 
S2 



(40) 



One defines an effective mass rries in the scalar sector: 

1 sin^ ip ^ cos^ ip 



m 



eflf 



One then has, for the K system, 



<2 



192t;2 m\ 



1 + 



m 



K 



12J 



11m 



K 



Hh 
12 



{nis + md 



12J 



1 + 



m 



K 



[ms + md) _ 



[(Nd);, + (Nd) 



12J 



1 + 



[m, + TTidY 



[{Nd);, - {Nd)J 



(41) 



(42) 



(43) 



Both rrix and fx are given in Appendix A. In equations (H2|) and (H3!) . we 
should note that the matrix Nd is real in our model, therefore both and 
M^'^ are real. 

In the K system, we use M12 to fit 



ArriK 



2\M- 



12 1 , 



Im(Mi2A^) 
^2 Amx |A„ 



2 ' 



(44) 
(45) 



where Xu = V*gVud- In the K system there are important long-distance contri- 
butions to M12, which we do not know how to compute precisely. Therefore, 
in that system we use for Mf^ only the short-distance box diagrams, but 
allow ArriK calculated by using equations (155]1 and (jSj) to be in between one 
half and twice the experimental value. 
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In the Bd and Bg system^ we fit Am^^ and Auib^ by using a formula 
analogous to equation 



There are uncertainties in the "bag parameters" used in Mf^. In Mf^, 
we use the vacuum- insert ion approximation to calculate the values of the 
hadronic matrix elements, and do not allow for corrections to the matrix 
elements provided by that approximation. In order to allow for these theo- 
retical uncertainties, we let our results for e^, Am^d, and AniBs differ from 
the experimental values by at most 10%. 

We fit two more quantities, sin (2/3) and sin {2a). These are computed in 
the following way. For the K, Ba, and Bg decays we defin^ 



q M, 



12 



P \Mu\ 

CP violation in B^ — )■ ipKs is determined by 



(46) 



By using equations (BLS-28.24), (BLS-30.34), and (BLS-30.35), we know 
that in the SM X^Ks = g^P (~2i/3), where /3 is a certain phase of the CKM 
matrix^ We therefore use 

sin(2/3) = -ImAv,i^^ (48) 

and compare our — Im A^i^g to the current experimental value of sin (2/3). In 
this way we constrain the NP contributions to M12 in both the Bd and K 
systems, through equations ( H6|) and ( P7|) . 

An isospin analysis of the decays B^ — vrvr may be used, together with 
the analysis of Bd — t- pvr and B^ pp., to extract 

K^^(^) (49) 

Using equations (BLS-28.24) and (BLS-30.35), we see that in the SM A,r7r = 
— exp {2ia). We thus use 

sin (2a) = — Im At^tt (50) 



®In those systems we use for Mf^ a simplified expression involving only the exchange 
of top quarks in the box diagram. 

^This definition uses the sign conventions in l^- Many authors use instead q — > —q. 
^°The phase e' in equation (BLS-28.24) is known to be tiny. 
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together with the current experimental value of sin (2a) to constrain in 
the Bd system. 

To summarize, we work with 14 parameters: a, b, c, x, y, a', b', c', x', 
y', 9, r = V2/V1, niA, and nies- With those 14 parameters we strive to fit 15 
observables: m„, rric, mt, rrid, rUs, rUb, \Vcb\, |Kfe|, Am^, ex, Am^^, 

AniB^, sin (2/3), and sin(2Q;). We found that the fit is possible and, indeed, 
we have found a large variety of input parameters, i. e. of points in parameter 
space, which are able to satisfy the criteria of the fit. 



4.2 Second stage 

Each one of the fits found in the previous subsection is a posteriori passed 
through a filter, to ensure that0 

the Yukawa couplings are perturbative; 

the quantity sin (2/3) computed from the decays Bd — D^D^ is correct; 
the angle 7 lies in the allowed range; 
Am£) is not too large. 

We next explain each of these four points. 

From equation f[T7|) we see that the Yukawa-coupling matrix Fi has matrix 
elements a/fi, and c/fi; likewise, the matrix F2 has elements x/v2 and 
y/v2. In order to preserve the perturbation expansion, we have required that, 
for any particular solution in our fit, all matrix elements of Fi and F2 — and, 
likewise, of Ai and A 2 — do not exceed 47r in modulus. 

In the decays Bd — )• D^D^ one ha^ 

?^ v,bv:d 



and 

sin(2/3) = ImAo+i?-. (52) 



^^We find that all the fit points which have passed through the filter actually have \Vtd\ 
in the SM range. 

Although there is a loop-suppressed (but not CKM-suppressed) penguin contribution 
to this decay, this can be ignored due to the large experimental error. It is only the sign 
of this observable which will be of use below. 
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We require sin (2/3) computed in this way to agree with the experimental 
value. Notice that this path to sin (2/3) does not include M12 in the K 
system. 

The CP-violating phase 7 = arg {—VudVchV*^y*^ has been extracted from 
the decays — )■ DK"^. There are two experimentally allowed regions: one 
region in which 7 70° is in the first quadrant and has values consistent 
with the SM, and another region with 7 70° — 180° in the third quadrant. 
The solutions with 7 in the third quadrant, though, are excluded by current 
measurements of the semileptonic asymmetry in B decays p]. We have 
computed 7 in each of our fit points and used it a posteriori in our fit. 

The experimental data discussed this far potentially constrain the scalar 
masses and m^.^ and the FCNC matrix N^- The most important con- 
straints on Nu come from mixing [i.e. M12) in the D system. In the SM, 
that mixing has three origins: box diagrams, dipenguin diagrams, and long- 
distance physics. The long-distance effects should be dominant but are very 
difficult to estimate reliably. Therefore, we only require that the NP con- 
tribution by itself alone should not exceed twice the experimental limit on 

We want to comment on a set of points that we have found at the first 
stage of our fit and which display an inverted unitarity triangle, i.e. have a 
negative Jarlskog invariant [10] Jckm = I^i (KtsKfeKTbKs)- Such points fit 
well the 15 observables used in the first stage, but are all eliminated at the 
second stage of the fit, because they display 7 ^ —70°, in contradiction with 
experiment. Besides, some of the Jckm < points suffer from the extra 
problem that they rely on dramatic contributions to M12 in the K system, 
with ~ —2Mf^. In these points the sign of q/p in the K system is 

inverted with respect to the SM. As a result, sin (2/3) extracted from ipKs 
decays would have the opposite sign to the sin (2/3) extracted from D^D~ 
decays, which is excluded by experiment. 

4.3 Two extra quantities 

CP violation has also been measured in the decay B^ — )■ ip(f). It is determined 



by 




(53) 
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Using equation (BLS-30.36), we see that the SM leads to A^,^ = — exp {2i(3s), 
where /3, is a phase in the CKM matrix which, in the SM, is of order a few 
percent o Thus, in the SM 

sin(2^,) = -ImA^,^. (54) 

We might have used the current measurement of sin (2/3^) from the decays 
Bs — )■ to constrain Mf^ in the Bg system. However, a recent average [11] 
excludes the SM at the 2.3 a level. Our fits always yield a Ps very close to 
its SM value; thus, our model does not provide a solution to this discrepancy 
between the SM and experiment. 

In this model, direct CP violation is negligible in D decays, and therefore 

m 

arg ( Tl, = 0, (55) 



relates to the amplitudes for the decays D — )• K^K^ . As a result, 

012 = arg (Mi2rt2) = - arg (^M*,^ ^f^^ . (56) 

As shown in reference [12], the theoretical parameter 0i2 can be extracted 
from the experimental data. 



5 The fit: results 

After the two stages of our fit we still have many points which have satisfied 
all the filtering criteria. With those points we have made a number of figures, 
which we next present. 

Figure [1] displays the asymmetry between and as a function of 
the smallest of those two masses. Clearly, if the scalar masses are both very 
large, then the model is effectively like the SM, except for the important 
fact that now the CKM CP-violating phase does not arise from complex 
hard (dimension-4) Yukawa couplings, as in the SM, but rather from a soft 
(dimension-2) CP-breaking term in the scalar sector. We find, however, that 
our model can have scalar masses as small as a few hundred GeV, especially 

^•^In [8] the phase (3s has been called e. 
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Figure 1: The masses of the scalars. We do not display the points where 
both rriA and rngfj are larger than 1 TeV. 

when rriA ~ ?7ieff0 In this case of low scalar masses, we have found that 
/Mf^ can be very large in the kaon sector, but is not larger than 10% 
in the Bd and Bg systems. 

In order to quantify the latter statement, we define [9] 

Mi2 = Mff + M[^/ = Mff A, (57) 

where the SM limit corresponds to A = 1. We shall use a subscript K,d,s 
in A to refer to the cases of the K system, B^ system, and Bg system, 
respectively. The current measurements do not agree well with the SM. 
Setting Ak = 1 and excluding the measurement of Ps, the CKMfitter Group 
|13j finds that the current constraints on A^ and A^ exclude the SM at 
the 2.2(7 and 1.9 a levels, respectively. The measurements of Ps are much 
above the SM prediction and further worsen this inconsistency [H]. Similar 
conclusions are drawn by the UTfit Collaboration [14]. 

Figures [2], [3l and H] contain the results of our fits for A^, A^, and A^, 
respectively. We see that Im A is in general quite small. This is a reflection 
of the fact that in our model CP violation lies exclusively in the CKM matrix 
while the matrix A*"^ is real; therefore Mf!^ is, in our model, real in all three 

^^Low scalar masses may in some cases be excluded by other experimental constraints 
that we have not taken into account, for instance by top-quark decays. 
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Figure 2: A parameter for the K mesons. 

neutral-meson systemslll We see in Figure [2] that Re Ak can be as large as 
three or four. This freedom is due to the large uncertainty in the long-distance 
contributions to K mixing. On the other hand, since ex is small, IirAk 
cannot be larger than two or three percent. In the Bd system, changes of 
of order 10% relative to the SM are possible both in the real and imaginary 
parts. For some of our points this decreases slightly the inconsistency of the 
SM with the experimental fits. However, this improvement is not dramatic 
because the experimental fits prefer ImA^ < and ReA^ < 1, while our 
points with Im A^ < have Re A^ > 1, cf. Figure [31 In the Bs system. Re A^ 
can differ from 1 by 10% or so, while Im A^ remains at the 0.1% level. Thus, 
in the Bg system our model is as (in) consistent with experiment as the SM. 

Figure [5] contains the predictions of our model for 0i2, based on the full 
set of our points and using exclusively M^^, i.e. assuming Mf^ = 0. We 
see in Figure |5] that sin^ 0i2 is, in our model, arbitrary; this illustrates how 
important CP violation in the D system can be in constraining models of new 
physics such as ours. Notice that the present experimental constraints 
on depend on a set of measurements which are highly correlated [13j : 

^^We have neglected potentially complex contributions to at loop level, notably 

box diagrams involving intermediate charged scalars C^. This is consistent with our 
previously stated assumption that the NP tree-level contributions to quark decays are 
much smaller than the SM ones. 
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Figure 3: A parameter for the mesons. 

precise numbers are not available, but we estimate, based on the method in 
[12], that sin^ 0i2 < 0.34 at the 1 a level. 

6 Conclusions 

In this paper we have presented a two-Higgs-doublet model with a sym- 
metry and the usual CP symmetry in all the hard (dimension-four) terms 
but broken in one, and only one, soft (dimension-two) term in the scalar po- 
tential. We have shown that this model displays a CP violation which, just 
as in the SM, is concentrated in the CKM matrix, even though it has a com- 
pletely different origin. Contrary to most other 2HDMs, our model exhibits 
CP violation neither in scalar-pseudoscalar mixing, not in the scalar self- 
interactions, nor in the matrices Nd^u which parametrize the flavour-changing 
Yukawa interactions of the neutral scalars. 

Our model has only eleven parameters — ten moduli and one phase — to 
fit the six quark masses and the four independent observables of the CKM 
matrix. When computing mixing in the neutral-meson-antimeson systems 
one needs three extra parameters — the ratio of VEVs, the mass of the pseu- 
doscalar, and a weighted mass of the two scalars. With these parameters one 
is able to fit most observables, just as in the SM. Remarkably, many of these 
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Figure 4: A parameter for the Bg mesons. 

fits display scalar masses as low as 400 GeV. 

We have emphasized the relevance that a measurement of CP violation in 
D-meson-antimeson mixing may have in eliminating some of our fit points 
and, thus, in reducing the viable parameter space of our model. 

Acknowledgements: The work of L.L. and of J.P.S. is funded by FCT 
through the projects CERN/FP/109305/2009 and U777-Plurianual, and by 
the EU RTN project Marie Curie: MRTN-CT-2006-035505. J.P.S. is grateful 
to Y. Grossman, Y. Nir, M. Papucci, and D. Pirjol for exchanges concerning 
solutions with negative Jckm- 



References 

[1] T.D. Lee, A theory of spontaneous T violation, Phys. Rev. D8 (1973) 
1226. 

[2] Papers published in the last few months include: 

P.M. Ferreira, H.E. Haber, M. Maniatis, O. Nachtmann, and J. P. Silva, 
Geometric picture of generalized-CP and Higgs-family transformations in 
the two-Higgs-doublet model. larXiv: 101(10935 [hep-ph]; 



18 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

min(m^, m^^) (GeV) 
Figure 5: Predictions for sin^</)i2. 

B. Grzadkowski, M. Maniatis, and J. Wudka, Note on custodial symmetry 
in the two-Higgs-doublet model, arXiv:1011.5228 [hep-ph]; 
H.E. Haber and D. O'Neil, Basis-independent methods for the two-Higgs- 
doublet model III: The CP-conserving limit, custodial symmetry, and the 
oblique parameters S , T, U, arXiv:1011.6188 [hep-ph]; 
B. Grzadkowski and P. Osland Natural two-Higgs-doublet model, arXiv: 
1012.0707 [hep-ph]; 

M. Jung, A. Pich, and P. Tuzon, The B — )■ rate and CP asymmetry 
within the aligned two-Higgs-doublet model, [arXiv:1011.5154l [hep-ph]: 
B. Grzadkowski, O.M. Ogreid, P. Osland, A. Pukhov, and M. Purmoham- 
madi. Exploring the CP-violating inert-doublet model, a rXiv: 1012.46801 
[hep-ph] . 

[3] P.M. Ferreira and J. P. Silva, Abelian symmetries in the two-Higgs-doublet 
model with fermions, larXiv: 1012.2874. [hep-ph] . 

[4] R.D. Peccei and H.R. Quinn, CP conservation in the presence of instan- 
tons, Phys. Rev. Lett. 38 (1977) 1440. 

[5] G. C. Branco and M. N. Rebelo, The Higgs mass in a model with two 
scalar doublets and spontaneous CP violation, Phys. Lett. B 160 (1985) 
117. 



19 



[6] P. M. Ferreira and J. P. Silva, A two-Higgs-doublet model with remarkable 
CP properties, Eur. Phys. J. C 69 (2010) 45. 

[7] K. Nakamura et al. (Particle Data Group), Review of particle physics, J. 
Phys. G 37 (2010) 075021. 

[8] G.C. Branco, L. Lavoura, and J. P. Silva, CP Violation (Oxford Univer- 
sity Press, Oxford, 1999). 

[9] A. Lenz, U. Nierste, and the CKMfitter Group (J. Charles et al.), 
Anatomy of new physics in B-B mixing, larXiv: 1008. 15931 

[10] C. Jarlskog, Commutator of the quark mass matrices in the standard 
electroweak model and a measure of maximal CP violation, Phys. Rev. 
Lett. 55 (1985) 1039. 

[11] D. Asner et al. (Heavy Flavor Averaging Group), Averages of b-hadron, 
c-hadron, and tau-lepton properties, larXiv:1010.1589t online update at 
http: / /www.slac. stanford.edu/xorg/hfag 

[12] Y. Grossman, Y. Nir, and G. Perez, Testing new indirect CP violation 
Phys. Rev. Lett. 103 (2009) 071602. 

[13] J. Charles et al. (CKMfitter group), CP violation and the CKM matrix: 
Assessing the impact of the asymmetric B factories, Eur. Phys. J. C 41 
(2005) 1, and online at http: / /ckmfitter.in2p3.fr/plots_ICHEP10[ 

[14] M. Ciuchini et al. (UTfit Collaboration), 2000 CKM-triangle analysis: 
A critical review with updated experimental inputs and theoretical param- 
eters, JHEP 07 (2001) 013, and online at http://www.utfit.org, 

[15] See, for example, A. Petrov, CP violation in charm, ' arXiv:1101. 38221 

[16] G. Buchalla, A.J. Buras, and M.E. Lautenbacher, Weak decays beyond 
leading logarithms. Rev. Mod. Phys. 68 (1996) 1125. 

[17] H.-W. Lin et al.. Charm as a domain wall fermion in quenched lattice 
QCD, Phys. Rev. D 74 (2006) 114506. 

[18] W. Grimus, L. Lavoura, O.M. Ogreid, and P. Osland, A precision con- 
straint on multi-Higgs- doublet models, J. Phys. G 35 (2008) 075001; The 
oblique parameters in multi-Higgs-doublet models, Nucl. Phys. B 801 
(2008) 81. 



20 



[19] J. Erler and P. Langacker, in [7]. 

[20] LEP Working Group for Higgs Boson Searches, Phys. Lett. B565 (2003) 
61; 

idem, Eur. Phys. Jour. C47 (2006) 547. 

A Input parameters 

We have used in our fits Gp = 1.16639 x 10"^ GeV"^ and mw = 80.4 GeV. In 
the neutral-kaon system, we have used tuk = 497.614 MeV, fx = 155.5 MeV; 
for the QCD correction factors of equation (BLS-17.16) we have taken jTB] 
rji = 1.38, ri2 = 0.57, and 773 = 0.47. For the "bag parameter" we have used 
[13] Bk = 0.723. Our results do not depend crucially on these inputs — small 
variations thereof do not change our conclusions. 

In the Bd system, we have used tub^ = 5.2795 GeV, /^^ = 190 MeV, 
TjB^ = 0.55, and = 1.219. For the Bg mesons, [HIE], tub^ = 5.366 GeV, 
/b, = 228 MeV, riB, = 0.55, and Bb, = 1.280. In the D system, = 
232 MeV [H] and = 1.86483 GeV [7]. 

We next present two of our fit points: one with low masses rriA and rries 
and another one in which one of the masses is low and the other one much 
larger. 

• First point: a = 105 MeV, b = 15.2 MeV, c = 6 MeV, x = 4.1387 GeV, 
y = 24.2 MeV, a' = 169.6184 GeV, b' = -8.5 MeV, c' = 1.2823 GeV, 
x' = 6.8821 GeV, y' = -1.8 MeV, 9 = 5.6548 rad, r = 0.4277, = 
415.6327 GeV, and rUes = 411.0434 GeV. 

• Second point: a = 17.5 MeV, b = 180.7 MeV, c = -27.8 MeV, x = 
4.2504 GeV, y = 74.9 MeV, a' = 172.8953 GeV, b' = -29.6 MeV, c' = 
-2.7 MeV, x' = -600.3 MeV, y' = 1.2920 GeV, 9 = 3.6987 rad, r = 
0.8217, mA = 400.9344 GeV, and m^s = 2.6596688 TeV. 

With these inputs one obtains the following values for the observables: 

• First point: = 5.8 MeV, = 107.8 MeV, = 4.1388 GeV, 
m„ = 1.8 MeV, = 1.2813 GeV, rrit = 169.758 GeV, |Ks| = 0.2256, 
\Vcb\ = 0.0405, \Vub\ = 0.0037, \Vtd\ = 0.0087, Atjik = 2.93 x lO'^^eV, 
AniB^ = 3.55 X 10"^ eV, Ams, = 1.13 x lO^^eV, je^] = 2.227 x lO'^, 
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JcKM = 3.113 xlO-^ sin (2a) = 0.1714, sin (2/3) = 0.71280 sin (2/3,) = 
-0.0397, 7 = 71.96°, = 1.3242 - 0.0032 i, = 0.9996 + 7.2444 x 
10~S, A, = 1.0000 - 9.6084 x IQ-^ i. 

• Second point: = 6.0 MeV, m, = 81.5 MeV, = 4.2542 GeV, 
m„ = 2.7 MeV, = 1.2923 GeV, rrit = 172.8963 GeV, |Ks| = 0.2249, 
\Vcb\ = 0.0425, \Vub\ = 0.0037, \Vtd\ = 0.0089, Am^ = 4.98 x lO'^eV, 
AmB, = 3.57 X 10"^ eV, Am^^ = 1.16 x lO'^eV, lex] = 2.128 x lO'^, 
JcKM = 3.211 xlO-^ sin (2a) = 0.1167, sin (2/3) = 0.74500 sin (2/3,) = 
-0.0407, 7 = 69.11°, Ak = 2.2183 - 0.0150 i, A^ = 0.9311 + 0.0678 
A, = 0.9088 - 0.0031 i. 



B Oblique parameters 

Relevant contraints on the scalar spectrum of a two-Higgs-doublet model 
arise from consideration of the so-called 'oblique parameters', especially of the 
parameters S and T0 Formulae for those parameters in a general MHDM 
have been presented in ref. [18]. In our particular 2HDM, one has 

T = — — —[cos'^i:f{ml,ml) + sm^ijf{ml;,ml) + f{ml,m\) 
— sin^ -ip f{ml, m\) — cos^ f{m\, m\) 

+ sin^ ^ f\ml) + cos^ ^ /'(m^) - f (m^)] , (Bl) 

where 

{X + y xy X 
In — X y, 
2 x-y y (B2) 
^ x = y, 

f'{m') = 3 [/(m|, m^) - /(m^^, m^)] . (B3) 

In equations fIBip and flB3p . mc is the mass of the charged scalars C^, mw 
is the mass of the W^, mz is the mass of the Z^, tjih is the mass of the SM 

^^This is the sin (2/3) which is obtained from the decays Bd — )■ ipKs- The value obtained 
from Bd D+D- is 0.7189. 

I'^The sin (2/3) obtained from Bd D+D- is 0.7487. 

^^The other obhque parameters are usuahy very small and, therefore, irrelevant. We 
have checked this explicitly for some of our points. 
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;iggs particle, and = 1 — m^/m|. The expression for S is 



4 2 



+ sin ?/'^(xi) + cos ^(x2) — ^'(a;^:/) + In ■ 



(B4) 



where Xk = mj^/m'^ for = 1,2, A, C. The functions g{x,y) and ^(a;) in 
equation ( 1B4p are in the second paper of ref. [I8jl^ 

For each of our fit points we only have the masses meg — in eq. pijl — and 
rriA- Equation (141 p may be solved for the mixing angle ip, yielding 

sin2^ = ^^ ^, cosV = ^^5 ^. B5 

m^g — mi m^g — m2 

Therefore, either mi < meg < m2 or m2 < m^s < mi. 

In order to check whether each of our fit points — defined by given values 
of mefi and — is compatible with the experimental bounds on the oblique 
parameters [19], we have inputted various values of mi^2 and mc- From mi 
and m2 we have computed ip through equation ( 1B5I) and then the oblique 
parameters S and T. With a fast fitting program we have been able to find, 
for a large part of our fit points, values of mi, mo, and mc such that 5" and T 
result compatible with the experimental bounds]^"! The masses mi 2,c can be 
chosen such that the parameter T does not result too large (either positive 
or negative). The parameter S usually turns out to be positive and relatively 
large {S > 0.1), but for mostly of our points it can be made compatible with 
the experimental bounds. 

As an example, one of our fit points has meir = 5.4711 TeV and m^ = 
679.9875 GeV. Choosing mi = 0.99 meg, m2 = 2.0283 mefi, and mc = 2 meg, 
one obtains 5* = 0.21 and T = 0.19. 



-•^^One has g{x,y) = G{xz,yz, z) and g{x) = G{xz,z), with the functions G{I,J,Q) in 
equation (C2) and G{I,Q) in equation (C5) of ref. [T8] . 

^°We have used niH = 117 GeV in accordance with one of the experimental elhpses in 
Figure 10.4 of ref [19]. 

^^We have not been able to explicitly find out, for all of our fit points, values of mi^2,c 
such that both S and T agree with the experimental bounds, but we cannot exclude that 
that is possible. 
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C Direct LEP bounds 



In Figured] we have shown that our fit sometimes yields scalar masses as low 
as 100 GeV. Since the current limit from LEP is 114.4 GeV [7], it is necessary 
to verify that our results do not contradict that bound. The LEP result is 
obtained by looking at the associated production of a scalar particle and a Z° 
boson, e^e~ — )■ Zh, which is possible due to the vertex ZZh. For a 2HDM 
there is also ZH production, but not ZA production, since there is no ZZA 
vertex. Moreover, as compared to the SM, the coupling of the vertex ZZh 
[ZZH) is reduced by factors related to the mixing angle if) in equation (l^Tj) . 
Indeed, in our model one has 



where g'zzs ~ sin^'?/' {cos^ tp) if S = h {S = H). Therefore, in a 2HDM it 
is possible to have scalars with masses lower than the LEP bound, provided 
those scalars couple more weakly to ZZ than in the SM. 

As explained before, our fit to the quark masses, CKM matrix elements, 
and CP- violating quantities has produced a large number of acceptable points 
in parameter space. Out of those, as seen in Appendix B, the vast majority 
conforms to the existing constraints on the oblique parameters. In Figure 
we plot, for h and H simultaneously, the comparison between the set of points 
which have passed the oblique-parameter fit and the experimental data from 
the direct searches at LEP; acceptable points must be below and to the right 
of the solid line in the plot. We see that, with the exception of only three 
points, the parameter space that we have found agrees perfectly with the 
LEP data. (As with the fit to the oblique parameters, we cannot exclude hat 
other values of mi 2,c can be found, such that all the points agree with the 
LEP experimental bounds.) 

We have also looked at the existing LEP bounds on scalar-pseudoscalar 
production. Those bounds extend to 225 GeV in the sum of the masses of the 
scalar and the pseudoscalar. We have found that all our points which survive 
the LEP bounds on Z°-scalar production display a sum of the masses of 
the scalar and the pseudoscalar which exceeds 225 GeV. Therefore, all those 
points also survive the LEP bounds on scalar-pseudocalar production. 
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Figure 6: Constraints on scalar masses from direct searches at LEP (data 
taken from ref. |20j). 
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